Metastable states influence on the magnetic behavior of the triangular lattice: 
Application to the spin-chain compound Ca3Co206 
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It is known that the spin-chain compound Ca3Co2 06 exhibits very interesting plateaus in the 
magnetization as a function of the magnetic field at low temperatures. The origin of them is still 
controversial. In this paper we study the thermal behavior of this compound with a single-flip Monte 
Carlo simulation on a triangular lattice and demonstrate the decisive influence of metastable states 
in the splitting of the ferrimagnetic 1/3 plateau below 10 K. We consider the [Co2 06]n chains as giant 
magnetic moments described by large Ising spins on planar clusters with open boundary conditions. 
With this simple frozen-moment model we obtain stepped magnetization curves which agree quite 
well with the experimental results for different sweeping rates. We describe particularly the out-of- 
equilibrium states that split the low-temperature 1/3 plateau into three steps. They relax thermally 
to the 1/3 plateau, which has long-range order at the equilibrium. Such states are further analyzed 
with snapshots unveiling a domain-wall structure that is responsible for the observed behavior of 
the 1/3 plateau. A comparison is also given of the exact results in small triangular clusters with 
our Monte Carlo results, providing further support for our thermal description of this compound. 

PACS numbers: 75.25. +z, 75.30.Kz, 75.40.Mg, 75.60.-d 



I. INTRODUCTION 

Low-dimensional interacting spins always reveal very 
interesting magnetic properties, as well as new electronic 
transport phenomena [K [2] [3] IJ. Along this line, sys- 
tems with geometric frustration have long attracted our 
attention because the ground- state properties, like de- 
generacy, are usually responsible for peculiar behaviors 
at low temperatures. Included are exotic magnetization 
dependences with external fields in hysteresis curves, for 
example. One group of compounds that exhibits such 
properties is the family CsCoXs, where X stands for CI 
or Br [5 . These materials are unidimensional Ising-like 
magnetic systems. The spin chains are arranged on a tri- 
angular lattice, while the intra- and interchain exchange 
couplings are both antiferromagnetic. 

Another fascinating spin-chain family of compounds 
with triangular arrangements has the formula A'aABOe 
(where A' may be Ca or Sr, whereas A and B are tran- 
sition metals). One particular system that caught our 
attention is CasCosOe E El [lOl [III [H [El [Til [IS 
[MIIIlliaillEQlEIlElESE^ which has a rombohe- 
dral structure composed of large chains of [Co206]n along 
the c-axis of a corresponding hexagonal lattice. The Ca 
atoms are located among those chains. The chains are 
made of alternating face-sharing CoOe trigonal prisms 
and CoOe octahedra. Each Co chain is thus surrounded 
by six equally spaced Co chains forming an hexagonal 
lattice in the a6-plane (see Figure [T]). 

Generally, the long Co-Co interchain distance (5.3 A), 



in Ca3Co206, compared to the short Co-Co intrachain 
distance (2.6 A), ensures the hierarchy of the magnetic 
exchange energies, Jintra ^ Winter, in modulus. This fact, 
together with a strong spin-orbit coupling [25l [26] , pro- 
vide a high uniaxial anisotropy along the c-axis. Neu- 
tron diffraction studies [7 , among others, have revealed 
a ferromagnetic intrachain interaction (Jintra < 0), while 
an antiferromagnetic interchain interaction (Jinter > 0) 
was observed in the a6-plane. Consequently, it is very 
plausible to describe the magnetism of this compound by 
considering an Ising triangular lattice of large magnetic 
frozen-moments linked to the chains, each one playing 
the role of a giant spin, pointing either up or down. 




FIG. 1: Hexagonal lattice of Co chains projected onto the ab- 
plane (blue points) of the spin-chain compound Ca3Co2 06. 
Each chain is ordered ferromagnetically and is coupled anti- 
ferromagnetically to its neighbors. 
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One of the interesting features observed in Ca3Co206 
is the appearance of plateaus in the magnetization curve 
M vs below Tq = 25 K, when the magnetic field is ap- 
plied along the chains (c-axis). For temperatures in the 
range, lOK < T < Tc, it is observed a rapidly increase 
of M with H. Close to zero field, the magnetization 
reaches a plateau at the value of Ms/ 3 (Mg, the satura- 
tion magnetization is ~ 4.8 /iB/f-u.). This value remains 
constant up to He ~ 3.6 T, where the magnetization 
springs up to its saturation value Mg. Otherwise, when 
T < 10 K, called the ferrimagnetic region, the observed 
plateau at M^/S, splits into three small steps, changing 
at fields H = 1.2 T, 2.4 T, 3.6 T, respectively. Much 
effort has been done to understand such behavior, and it 
is still a controversial matter. Some papers [161 1211 1211 
roused the question that this behavior might be related 
to experiments involving single-molecule magnets (SMM) 
[27| 28 , where a Quantum Tunneling of the Magnetiza- 
tion (QTM) is realized [29 . We restrain such statement, 
at least for temperatures above 4K, as we demonstrate 
here and leave the question of QTM effects to be searched 
below 2K, according to the experimental situation. 

We follow instead the line of reasoning that metastable 
states, or peculiar mozaic configurations, might play a 
role in the splitting of the magnetization of this frustrated 
material, as first suggested by Maignan et al. [9 . There- 
fore, we started a simulation along the lines given by Ku- 
dasov [30 , who predicted an opening of the 1/3 plateau 
using Monte Carlo analysis at T = 0. More recently, 
Yao et al. ^TJ |32] , studied this compound with numer- 
ical Monte Carlo simulations on finite clusters. The lat- 
ter works describe the system with an Ising Hamiltonian 
on a triangular lattice, using periodic boundary condi- 
tions. Despite their results reproduce the steps width, 
AH ~ 1.2 T, in the magnetization curves, they do not 
find a clear-cut {equilibrium) configuration that explains 
the splitting of the 1/3 plateau. In this work we report 
on a similar study using the same model, but consider- 
ing instead free boundary conditions on an anisotropic 
system, which is added for completion with the intrinsic 
dipolar term at the end. We apply non-equilibrium tech- 
niques on this system. We detect in this form, by sweep- 
ing the field at different rates, the presence of metastable 
domains- walls that perfectly explains the observed trends 
of the splitting of the 1/3 plateau at low temperatures, 
and verify that the latter has long-range order at the 
equilibrium. To be consistent with the previous exist- 
ing models, we use the same set of parameters for the 
Hamiltonian [^0^, 211 El] , namely, an interchain antiferro- 
magnetic interaction, J = Jinter = 2.25 /ieV, and a large 
Ising magnetic moment, S = 32. Disorder or dispersive 
effects in the exchange coupling constant were not re- 
quired to obtain the splitting of the 1/3 plateau, as it 
will become clear in the following. Therefore, our model 
and methods are different from Yao et al. [3TJ |32 , and 
we concentrate our efforts mainly to describe the effects 
of the sweeping rate on the observed metastable states, 
that generate the sub-steps in the magnetization. 



II. MONTE-CARLO SIMULATIONS IN 2D 

The spin-chain compound Ca3Co206 is formed by 
close-packed one-dimensional chains along the c-axis, the 
chains span a triangular lattice over the a6-plane. Due 
to such geometric configuration and the hierarchy of the 
exchange interactions, the system can be modeled as a 
two-dimensional (2D) antiferromagnetic Ising model with 
nearest-neighbor exchange interactions [30l |31] |32] . The 
magnetic moments are representative of the [Co2 06]n 
chains, which are always perpendicular to the a6-plane. 
The choice of this frozen-moment model is a consequence 
of a time-scale assumption. We believe that the charac- 
teristic time involved in reversing the individual spins of 
each chain is much shorter than the rearrangement of the 
triangular lattice by changing the magnetic field. This 
hypothesis is based on previous numerical results [33 on 
nano-wires, as well as on the large number of sites that 
are involved in changing the whole triangular lattice as 
compared to those on a single spin-chain. 

We thus performed Monte Carlo simulations using a 
triangular Ising model, where the magnetic moments are 
frozen, pointing either up or down, and coupled antifer- 
romagnetically (J > 0) to their nearest neighbors. In 
this picture, the Hamiltonian H of the system is 

H = J^S^S]-g^iB^S,-U, (1) 

where is the spin of each magnetic moment, H is the 
external magnetic field along the +c-axis, ^ = 2 is the 
electronic Lande factor, jub the Bohr magneton and (z, j) 
indicates summation over all pairs of nearest-neighbors 
on a triangular lattice. We used a single-fiip Monte 
Carlo method together with the standard Metropolis al- 
gorithm to reject/accept fiips in the magnetization with 
a Boltzmann factor [34]. Magnetization measurement 
at any given field was done using 1,000 averages over 
a fixed number of accepted Monte-Carlo (MC) fiips. We 
used different 2D finite cluster configurations, including 
super- hexagons, octagons, dodecagons and other more 
circularly-shaped geometries to represent the samples. 
The magnetization curves found in all these finite clusters 
were almost the same, the main difference was the num- 
ber of MC fiips needed to reach the plateaus for different 
sizes. Due to space limitations and visualization aspects, 
a selection of the results for one cluster size with hexag- 
onal geometry was chosen. Extra material showing the 
evolution of magnetic configurations with the field can 
be found at the website by the interested reader [supple- 
mentary material]. 



III. MAGNETIZATION RESULTS 

In this section we discuss our magnetization results 
by using Monte Carlo calculations on finite clusters with 
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FIG. 2: MC results for the magnetization as a function of the 
appUed field for three different temperatures in the ferrimag- 
netic phase, T < 10 K, and N = 2431 sites. The splitting 
of the 1/3 plateau into three successive plateaus is observed 
below 10 K. They extend up to fields H ^ 1.2 T, 2.4 T and 
3.6 T, respectively. We see that they merge into a unique 
plateau at M ^ Ms/ 3 as T increases (see arrows). 



FIG. 3: MC results for the magnetization as a function of 
the applied field for three different sweep rates, at T = 4 K, 
and N = 2431 sites. This is accomplished by letting more or 
less accepted MC flips (mcf) per magnetic field. Notice the 
tendency of merging into the plateau Ms/ 3 as the number 
mcf increases (see arrows), indicating a time decay towards 
the 1/3 state. This behavior will be further analyzed. 



open boundary conditions, always increasing the field, 
starting at = T. The initial magnetic state was ob- 
tained by thermalizing a random disordered configura- 
tion, performing 10^ MC flips in a training period. After- 
wards we initiated the increase of the field at a constant 
sweep rate. We calculated clusters with hexagonal ge- 
ometries for N = 1387, 2431, 3169, 4921 and 6931 sites. 
Some other larger clusters were partially calculated, but 
we verified that the agreement with the experimental re- 
sults did not improve with size. As the system size in- 
creased the plateaus structure remained unaltered, only 
small perturbations (sub-steps) on each of the three main 
plateaus tend to disappear. 

In Figure [2] we readily observe the formation of three 
well defined plateaus below T = 10 K. Important to see is 
that they merge into one, at Ms/3 as T increases, indicat- 
ing a close resemblance with the experimental situation 
(see e.g. Fig. 3 in [E]). The presence of extra tiny sub- 
steps is an unavoidable boundary effect due to the finite- 
ness of the cluster used (in this case, N = 2431 sites). 
They are created by border isolated atoms whose con- 
tribution is strongly suppressed for larger cluster sizes, 
therefore, irrelevant for our argument. Another impor- 
tant feature to note in the results of Figure |2j shared 
by experiments, is the steeper rise of the magnetization 
for lower temperatures at the end of the 1/3 plateau 
{H,^3.6 T). 

Now, as a way to check how stable are such plateaus, 
we have performed MC analyses with an increasing num- 
ber of steps per magnetic field. In fact, as the MC sim- 
ulation is done with the accepted configurations, what 
really counts is MC flips (mcf) and not MC steps (mcs). 
We clearly see, in Figure |3j the tendency to merge the 



three plateaus into one as a function of the number of 
mcf. This certainly must be related to the experimen- 
tal sweeping rate [19 . This connexion allows us to say 
that the three plateaus are metastable states and are, 
therefore, a dynamical effect; for this very reason the ex- 
perimental magnetization curves are strongly dependent 
on the sweeping rate [19]. 

This merging effect can be further explained in the 
graph of Figure |4j where we have produced an almost 
perfect superposition between magnetization curves for 
different temperatures and, at the same time, different 
mcf. In such a case, we tried to draw a parallel between 
the mcf in our simulation and the sweeping rate in the 
experimental results, where this behavior is observed. It 
is obvious from this result the strong interplay between 
the effects of time and temperature (see e.g. Fig. 4 in 
[19]), a feature that was previously used as a fingerprint 
of QTM. Although mcf are not related to lab time mea- 
surement, we can attribute them to some sweeping rate 
by comparing the formation of the plateaus. Hence, the 
three sub-plateaus are made of metastable states, which 
develop in time to the 1/3 plateau, as in Figure [sj 

Dipolar interactions were also considered, as an in- 
trinsic effect always present in nanowires and other 
low-dimensional systems, but its overall effect here was 
merely to induce an effective antiferromagnetic field that 
retarded the field at which the steps occur. The dipolar 
term T^dip that was added to equation ([T]) is given by: 

nj 1 V- mj nij -3(m^-r^^-)(m^--r^^-) 

^dip = 2 2^ ^. 

where the summation is over all sites, and = i/iB^i. 
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FIG. 4: Magnetization versus applied field for two different 
temperatures and different sweep rates. We see that they 
merge into very similar curves. This overlap indicates that 
they correspond to the same statistical configuration, one 
which realizes similar magnetic moments. Such scaling is 
an extra indication of metastable configurations. Differences 
are seen at both steep rising regions, close to zero and at 
He ~ 3.6 T, fact which is also seen in the experiments ^19j . 

We see such effect in Figure [5| with the dipolar inter- 
action playing the role of an effective antiferromagnetic 
field, H ^ H -\- H^^ added to equation ([T]), that moves 
the magnetization curves to the right. We noticed that 
adjusting the exchange constant to 0.9 J, we can almost 
reproduce the results obtained from our previous 2D Ising 
model, without dipolar interaction (see Figure |5|. For 
this reason, we think, the magnetization curves are well 
described using only the exchange and Zeeman terms. 



IV. SNAPSHOTS AND DOMAIN- WALLS 

The physics behind these results can be visualized by 
using raw snapshots of selected configurations. In the 
next pictures we show different snapshots of the magneti- 
zation curves depicted with the red line of Figure [3] They 
were taken at three different magnetic fields: H = 0.81 T, 
H = 2.0 T, and H = 2.98 T, representative of each of the 
three sub-steps, into which the 1/3 plateau is splitted at 
low temperatures (T < 10 K). 

Our color convention to analyze the snapshots of the 
subsequent examples is depicted in Figure [6j Following 
[9^, we define four triangles and color them differently 
depending on the magnetization that results by summing 
the three spins at their vertices: see Figure |6] 

Configurations at zero field are formed by disordered 
bicontinuous patterns of small compact regions of violet 
and green phases, with total zero magnetization. For a 
field close to 0.25 T, domain walls, as seen in Figure [7| 
stabilize their motion to quasi-static domain walls until 
the field reach values of about 1.2 T. 
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FIG. 5: Magnetization versus applied field considering dipolar 
interactions, at T = 4 K, and N — 2431 sites. The red 
curve is for the original J value used in Figure |3] without 
dipolar interaction. The blue curve is for the same J, but 
now with dipolar interaction, and the green curve is for a 
reduced exchange, 90% of the original value, and the same 
dipolar interaction of the blue curve. The scaling of the red 
and the green curves is not perfect, but very close. 



Moving the field across the boundary, i7 = 1.2 T, pro- 
duces a rapid rearrangement of the linked domain walls 
positions, by changing simultaneously their internal con- 
figuration. They are now formed by mixed violet and red 
triangles, as shown in Figure [8] They rate together to 
zero extra magnetic moment, sometimes oscilating above 
and below, depending slightly on the simulation param- 
eters. The magnetic configuration in Figure [S] averages 
to zero, leading to M = M5/3, the value of the green 
ordered zone (formed by the hexagonal unit cells). 

Finally, for fields beyond the other boundary, H = 
2.4 T, the situation is quite different. Now the linked 
domain walls are formed mostly by red colored triangles 
(Figure 9) with a tendency to coalesce in small groups, 
or seeds, at the borders of the simulated clusters. This 
interesting feature of red seeds at the boundaries, as seen 




FIG. 6: The magnetic moment of each triangle depends on the 
orientation of each ferromagnetic chain. The four possibilities 
are schematized by different colors, but M = — 3 is seldom 
seen for positive magnetic fields in the following snapshots. 



FIG. 7: Snapshot of a particular magnetic configuration in an 
hexagonal cluster of 1387 spins, at the first sub-step, of field 
H = 0.81 T, at T = 4 K and 10^ mcf, as depicted with a 
red line in Figure [3] Notice the ordered substrate green zone, 
with M = Ms/3, superimposed by interlinked domain walls, 
formed mostly by connected triangles in violet color. These 
domain walls contribute to a reduction of the Ms/3 value of 
the green zone to M 0.22Ms, in this example. 



FIG. 8: Same parameters as in Figure [7| except at a higher 
field H — 2.0 T, corresponding to a magnetization value of 
0.33Ms in this case. Observe now the green extended zone 
superimposed by linked domain walls of mixed violet and red 
colored triangles, in a proportion that contribute almost zero 
extra moment. These linked domain walls, as those before, 
percolate throughout the system and are rather stable across 
a finite window of the magnetic field: 1.2T<i7<2.4T. 



V. TRIANGULAR CLUSTERS 



in Figure |9| is a peculiar phenomena very similar to pre- 
nucleation phases, as found in irreversible statistical me- 
chanics studies. Proper nucleation really starts beyond 
the critical field, He = 3.6 T, in this way: first coales- 
cence of small red seeds at the borders and then myriads 
of red nuclei forming and growing, taking account of the 
whole cluster, in the bulk, until saturation. The nucle- 
ation process above He is a very rapid, disordered one, 
destroying the linked domain walls structures seen in the 
former examples. Perfect alignment is only achieved at 
the end, for fields higher than 4 T. 



All these non-equilibrium phenomena just described 
were obtained for a particular sweep rate, of 10^ mcf. We 
tested various sweep rates, observing similar behavior for 
many of them. But, none of these metastable interlinked 
domain walls were present, however, when we go beyond 
a sweep rate of 10^ mcf. In that case, the disordered 
initial phase goes immediately to the 1/3 plateau, and 
beyond He = 3.6 T, to saturation. Such 1/3 plateau is 
represented solely by the green zone, without any further 
domain wall, thus of long-range nature. Therefore, we 
believe our characterization of such metastable states, 
sweep-rate dependent, must be related to the multi-step 
phases of Ca3Co206, found experimentally at low T. 



In this section, we show how the 1/3 plateau of the 
Ising model of previous sections shares interesting simi- 
larities with the simplest small clusters of triangular ge- 
ometry. A generic spin Hamiltonian for these clusters is 
the following Heisenberg Hamiltonian: 
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H = i ^ J., S, • S, - (5 ^ (S, • n) ' - ^/iB £ S, • H (3) 



Results for a triangular molecule are shown in Figure [To] 
M (the thermal average of the total spin z-component) is 
zero only at zero field, and it increases rapidly to the 1/3 
plateau for any finite H, crossing through H = hi until 
H = where it rapidly increases again to saturation. 
Such a transition is first-order at T = only. Increasing 
the temperature produces a smooth growing up of M 
due to thermal fluctuations. Notice that larger thermal 
fluctuations accelerate this trend above hi and retard it 



above as seen from Figure 10 



Interesting to note is that such a trend in Figure [TO 
resembles the MC solution of the Ising model in the tri- 
angular lattice of our previous sections. Although it is 
not the same behavior, a tendency to form a 1/3 plateau 
is without any doubt. No further sub-plateaus can possi- 
bly be observed in a triangular molecule because of lack 
of phase space. Domain walls in the Ising dynamics are 
minimally formed by triangles, as seen before. 



FIG. 9: Snapshot of the third sub-plateau, at a field H = 
2.98 T. Other parameters as in Figures|7|and[8] The magnetic 
configuration is now formed by the green extended region of 
M = Ms/3 superimposed by domain walls in red color, giving 
a total magnetization of 0.41Ms, in this case. Notice also 
some red seeds distributed at the borders of the cluster, like 
in pre-nucleation phases of irreversible studies. Nucleation 
really starts beyond the critical field. He = 3.6 T. 
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FIG. 10: Exact magnetization against a magnetic field for 
the Heisenberg triangular molecule (A^ = 3), and for different 
temperatures, = feT, as shown in the legend. We see a 
tendency to form a 1/3 plateau. Other parameters used are: 
J = 2.25 /ieV, 5 = 1 /xeV, and 5' = 1/2. 



Further supporting arguments come from numerical 
exact solutions for (larger) finite clusters. In Figure 
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see, for example, numerical results for the magnetization 
of the pure Heisenberg model in the triangular lattice, in 
the absence of anisotropy, at T = 0, for = 36 and 39 
sites [38]. The formation of a clear 1/3 plateau of finite 
width, Ai7, is seen from the graphic. Such results have 



FIG. 11: Exact magnetization against a magnetic field for 
Heisenberg triangular clusters at T = 0, ^ = 0, and S = 1/2. 
The green and red lines are for N — and 39 sites, while 
the bold blue line is an extrapolation to the thermodynamic 
limit (see text). Notice the 1/3 plateau of finite width /\H 



been taken in the literature [39^, '401 El HI] as the evi- 
dence of a definite plateau in the thermodynamic limit: 
the bold line in Figure [TT] was drawn by connecting the 
midpoints of the finite-size steps for A/" = 39 sites, except 
at the boundaries of the 1/3 plateau. A good estimate for 
the lower and upper fields of the boundaries of the exact 
1/3 plateau in the triangular lattice were obtained com- 
paring spin- wave and exact diagonalization results. They 
are given by = 1.378 J and Hu = 2.155 J, respectively. 
See more recent results in [42^. 

To better understand our point, we should recall that 
both models, Ising and Heisenberg, show a similar be- 
havior - with a 1/3 plateau - because the magnetic field 
already breaks the symmetry from SU{2) to t/(l), thus 
making the strong quantum fluctuations inoperant in the 
isotropic Heisenberg triangular case and, consequently, 
the M = 1/3 state becomes a classical, collinear state, 
with long-range order. 



VI. CONCLUSION AND FINAL REMARKS 

In summary, we did an extensive investigation of the 
triangular lattice using Monte Carlo simulation tech- 
niques. We studied this lattice using a 2D Ising model for 
large magnetic moments coupled antiferromagnetically. 
We were motivated by multiple steps observed in the 
magnetization of the spin-chain Ca3Co206 compound. 
We found a very good agreement between our results 
and the experimental situation [19 . We demonstrated 
that the observed plateaus of the ferrimagnetic phase, 
below T < 10 K, consist of met ast able states and for this 
reason are time-dependent. In the experimental results 
such dependence was observed in the sweeping rate of the 
magnetic field, and in our case this dependence comes 
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from the number of Monte Carlo flips (mcf). In fact, a 
close parallelism is found between these two approaches, 
as seen from our Figures |2j |3j |4j and Figures 3, 4, 5 from 
Ref. [19 . Besides such agreement, we also found that 
in the limit of very large number of mcf, namely, a very 
slow rate of the fleld, the three sub-steps converge into 
one plateau, at M/Mg = 1/3, like in the experiments, 
thus giving a thermal decay of such metastable states 
into the equilibrium, flnal 1/3 plateau, which has long- 
range hexagonal order (described by the green ordered 
zone of the snapshots, with an hexagon as unit cell). 

Further analysis using snapshots of some speciflc 
conflgurations revealed the presence of linked (mobile) 
domain- walls that give origin to the splitting of the 1/3 
plateau below T < 10 K, as explained in the text. Such a 
result disproves the possibility of quantum tunneling as 
generating the multiple steps, at least for temperatures 
above 4 K, and for the reported sweep rates. Based on 
our results, we believe that at lower temperatures an ar- 
rested conflguration may provide a slower dynamics of 
the domain- walls, that would explain a T-independent 
time-decay process below 4 K, by critical slowing down, 
an aspect that certainly must still be checked down. 

Our solution has been checked by further considering 
the intrinsic dipolar interactions, flnding no relevant con- 
tribution from such couplings. We therefore discarded 
them as possible contribution of hysteresis phenomena 
which are strongly present at the ferrimagnetic region, 
below 10 K, in the experiments [19 . We believe that 
other ingredients and methods, like 3D simulations, must 
be considered to deal with such irreversible aspects, but 
this is beyond our model of the magnetization plateaus. 

Another point that we wanted to emphasize was linked 
to exact solutions of the Heisenberg model in the trian- 
gular lattice for flnite clusters. We did a comparison with 
our previous numerical Ising solution and found support- 



ing arguments with the presence of the 1/3 plateau in 
the isotropic Heisenberg case at T = 0. Both models 
share the 1/3 plateau, although for different reasons. In 
the quantum-mechanical case the applied magnetic fleld 
strengthens the frustration effects in the triangular lattice 
up to a point of having sufficient overlap with the clas- 
sical, collinear, three- fold degenerate, 1/3 plateau [41]. 
The point here is that the presence of the 1/3 plateau in 
the isotropic quantum limit does not require a strong uni- 
axial anisotropy. Although we did not a thermal study of 
the exact solution of the Heisenberg model, a feature hard 
to be done, we think that at low temperatures the 1/3 
plateau would survive to thermal ffuctuations in order to 
be observed numerically. Such a comparison, therefore, 
strongly supports our thermal, classical study of the spin- 
chain compound, Ca3Co206. 

While this paper was considered for publication we be- 
came aware of another work, by Kudasov et al. [43 , on 
the same problem, namely, a study of the dynamics of 
magnetization in frustrated spin-chain system Ca3Co206 
using a 2D model. Although they present similar results 
for the magnetization steps, they are complementary to 
ours in some respects. 
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